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This paper presents a theoretical investigation on kinematic and thermodynamic analysis of a beta type 
Stirling engine with rhombic-drive mechanism. Variations in the hot and cold volumes of the engine were 
calculated using kinematic relations. Two different displacer cylinders were investigated: one of them 
had smooth inner surface and the other had axial slots grooved into the cylinder to increase the heat 
transfer area. The effects of the slots grooved into the displacer cylinder inner surface on the performance 
were calculated using nodal analysis in Fortran. The effects of working fluid mass on cyclic work were 
investigated using 200, 300 and 400 W/m 2 K convective heat transfer coefficients for smooth and grooved 
displacer cylinders. The variation of engine power with engine speed was obtained by using the same 
convective heat transfer coefficients and isothermal conditions. The convective heat transfer coefficient 
was predicted as 104 W/m 2 K using the experimental results measured from the prototype engine under 
atmospheric conditions. The variation in cyclic work determined by the experimental study was also 
compared with the theoretical results obtained for different convective heat transfer coefficients and iso¬ 
thermal conditions. 

© 2013 Elsevier Ltd. All rights reserved. 


1. Introduction 

The Stirling engine is an external combustion machine, in¬ 
vented by Robert Stirling in 1816. The thermodynamic cycle of 
the Stirling engine consists of two isothermal and two constant 
volume processes [1,3-5]. Since the Stirling engine is externally 
heated, it can be operated with solar radiation, geothermal energy, 
natural gas, biomass, fossil fuels and all kinds of energy sources 
[2,6,7]. In Stirling engines, mostly air, helium and hydrogen are 
used as the working fluid [8,9]. Stirling engines have some advan¬ 
tages such as high thermal efficiency, low noise, low pollutant lev¬ 
els, and long life when compared with internal combustion engines 
[ 8 - 11 ]. 

Kinematic Stirling engines are classified into three different 
mechanical arrangements as a, |3, and y [3,5,12]. oe-type engines 
use two pistons situated with a 90° phase angle in separate cylin¬ 
ders. Cylinders can be configured in the form of a V or placed as 
parallel to each other. In (3-type Stirling engines, the cycle is con¬ 
ducted by a piston and a displacer operating concentrically in the 
same cylinder. In y-type engines, a piston and a displacer are situ¬ 
ated in separate cylinders. The displacer cylinder provides heating 
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and cooling of the working fluid. A regenerator can be placed inside 
or outside the displacer [3,12-14]. 

Different driving mechanisms such as a crankshaft, Ross-yoke, 
lever-drive, rhombic-drive, Scotch-yoke or swash-plate are used 
in Stirling engines. The cycle obtained by these driving mecha¬ 
nisms is different from the ideal Stirling cycle depending on the 
mechanical configuration. Cyclic work is the most important 
parameter in engine design [15]. Thermodynamic analyses are 
quite significant in determining the engine performance and 
parameters that affect the performance. Various analysis methods 
have been developed since the invention of Stirling engines and 
engine performance has been estimated during the design process. 

The first analysis for determining the thermodynamic perfor¬ 
mance of Stirling engines was developed by Gustav Schmidt in 
1871 and thermodynamic analyses of a, |3, and y type engines were 
conducted [3,6]. The nodal analysis method was developed in 1967 
by Finkelstein. In this analysis, the total volume of the engine was 
divided into 13 nodal volumes: the expansion volume (1), heater 
volumes (3), regenerator volumes (5), cooler volumes (3), and cold 
volume (1). Instantaneous temperatures, system pressure, cyclic 
work, temperatures and thermal efficiency are calculated in the 
nodal volumes [15-17]. In 1978, Martini developed a five-zone no¬ 
dal-isothermal analysis. In this model, the temperature of the 
working fluid in the nodal volumes was assumed to be equal to 
the wall temperature [18]. In 1994, Ladas and Ibrahim conducted 
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Nomenclature 



At 

nodal values of heat transfer surface (m 2 ) 

A rrij 

variation of nodal mass within a time step (kg) 

c p 

specific heat at constant pressure (J kg -1 K -1 ) 

R 

gas constant (J kg -1 K -1 ) 

Cy 

specific heat at constant volume (J kg -1 I< -1 ) 

Rcr 

radius of crankshaft (m) 

Et 

enthalpy flow in and out of the nodal volumes (J) 

T, 

nodal values of working fluid temperature (K) 

e c 

length of cold volume (m) 

Tc 

cold-end temperature of displacer cylinder (I<) 

fh 

length of hot volume (m) 

Th 

hot-end temperature of displacer cylinder (K) 

h p 

9 

distance between piston top and piston rod top (m) 

T„j 

nodal values of heat transfer surface temperature (I<) 

h 

displacer length (m) 

t 

time (s) 

h, 

nodal values of convective heat transfer coefficient 

At 

period of time steps (s) 


(Wm 2 K ') 

Uc 

length from cylinder top to rhombic center (m) 

^dr 

length of displacer rod (m) 

Vi 

nodal values of volume (m 3 ) 

Ipr 

length of piston rod (m) 

fir 

angle made by connecting rod with vertical (rad) 

Irh 

length of displacer and piston connecting rods (m) 

6 

total crankshaft rotation (rad) 

m t 

total mass of working fluid (kg) 



m, 

nodal values of working fluid mass (kg) 




a finite-time thermodynamic analysis of a Stirling engine. They 
analyzed the effects of regeneration and engine speed on engine 
power and engine efficiency by using the conservation law of mass 
and the energy equations [19]. In 2000, Kaushik and Kumar deter¬ 
mined the thermodynamic performance of a Stirling engine using 
finite-time thermodynamics. They investigated the effects of the 
temperature of different heat exchangers and heat source on en¬ 
gine power and engine efficiency [20]. In a study conducted by Kar- 
abulut et al. (2006) the power cylinder was concentrically situated 
in the displacer cylinder to increase the specific power and de¬ 
crease the external volume of a gamma type Stirling engines. The 
thermodynamic performance characteristics of the engine were 
predicted by a nodal analysis program in Fortran [21]. Kongtragool 
and Wongwises (2006) investigated the performance of a Stirling 
engine with the dead volumes of the cold volume, hot volume 
and regenerator using an isothermal model [22]. Parlak et al. 
(2009) conducted a thermodynamic analysis of a gamma type Stir¬ 
ling engine using a quasi-steady flow model. The analysis was per¬ 
formed for five volumes: compression, expansion, cooler, heater 
and regenerator. The variations of temperature, pressure, mass 
and work were calculated by using the conservation law of mass 
and the energy equations [23]. Timoumi et al. (2008) investigated 
the effects of regenerator sizes and materials on the performance 
of a General Motors GPU-3 Stirling engine [24]. Eid (2009) con¬ 
ducted a thermodynamic analysis of a beta type Stirling engine 
with a regenerative displacer by using the Schmidt analysis meth¬ 
od [13]. Cheng and Yu (2010) carried out the thermodynamic anal¬ 
ysis of a beta type Stirling engine with a rhombic-drive mechanism 
by using a numerical model. The energy equations for control vol¬ 
umes were derived and solved by taking into account the non-iso- 
thermal effects, the effectiveness of the regenerative channel and 
the thermal resistance of the heating head [12]. 

In this study, kinematic and thermodynamic analysis of a beta 
type Stirling engine with a rhombic-drive mechanism was con¬ 
ducted using the nodal analysis method. In the analysis, the in¬ 
crease of the heat transfer surface area inside the displacer 
cylinder on engine performance was investigated. The perfor¬ 
mance of the engine was estimated for two different displacer cyl¬ 
inders, of which one has a smooth inner surface and the other has 
axial slots grooved into the cylinder. Engine performance parame¬ 
ters were also calculated for different values of the convective heat 
transfer coefficient, working fluid mass and engine speed. The the¬ 
oretical results were compared with experimental data obtained 
from the prototype engine. 


2. Test engine 

Due to the lower side thrust and more silent operation [6], the 
rhombic drive mechanism was preferred and the engine was de¬ 
signed as beta type. The cycle is conducted by a piston and a dis¬ 
placer operating in the same cylinder. The rhombic drive 
mechanism consisted of a piston yoke, a displacer yoke, two pis¬ 
tons and two displacer connecting rods with equal length hexagon 
and two symmetric spur gears of equal diameter (Fig. 1). The spur 
gears rotate in opposite directions. The displacer rod is connected 
to the displacer yoke at the lower end of the equal length hexagon; 
the power piston rod is connected to the power piston yoke at the 
upper end of the equal length hexagon. The displacer and power 
piston yokes have linear motion. A displacer and a power piston 
rods are connected to pins on spur gears. The cold volume is be¬ 
tween the power piston and the displacer, and the hot volume is 



Fig. 1 . Beta type Stirling engine with rhombic-drive mechanism. 
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Table 1 

Technical specifications of the engine. 


Engine type 

P 

Driving mechanism 

Rhombic 

Diameter of displacer cylinder (mm) 

86 

Diameter of displacer piston (mm) 

84.5 

Length of displacer (mm) 

155 

Piston diameter (mm) 

86 

Length of piston (mm) 

90 

Rhombus length (lrh) (mm) 

66 

Radius of crankshaft (mm) 

4/2.6666 

Displacer rod diameter (mm) 

16 

Hot source temperature (T H ) (K) 

773 

Cold source temperature (T c ) (K) 

300 

Working fluid 

Air 

Engine speed (rpm) 

500 


above the displacer piston. The gap between the displacer and its 
cylinder is used as the regenerator. The schematic illustration 
and technical specifications of the engine are given in Fig. 1 and 
Table 1. 

As seen in Fig. 1, where the spur gear pin is at A, working fluid is 
compressed in the cold volume and the hot volume is at its mini¬ 
mum. At this position of the spur gear pin 9 is assumed to be zero. 
While the spur gear pin moves from A to B, the power piston re¬ 
mains around the bottom dead center of its stroke. The displacer 
moves downwards and working fluid is heated up by being 
pumped from the cold volume to the hot volume. The reduction 
in cold volume and expansion of hot volume are equal to each 
other, so the total volume remains constant. Thus, this process is 
carried out at a constant temperature. When the spur gear pin 
moves from B to C, the power piston and displacer move simulta¬ 
neously downward. While the cold volume remains constant at a 
minimum value, the working fluid expands in the hot volume by 
producing work. In this process, a lot of working fluid expands in 
the hot volume with the use of the rhombic-drive mechanism. 
When the spur gear pin moves from C to D, the power piston con¬ 
ducts a small displacement around the bottom dead center of its 
stroke. The displacer moves to the top dead center of its stroke. 
The displacer expels the working fluid from the hot volume to 
the cold volume, performing a cooling process at constant volume. 
When the spur gear pin moves from D to A, the displacer moves to 
the top dead center of its stroke and the hot volume is reduced. The 
power piston moves from down to up and the working fluid is 
compressed in the cold volume. Thus, the cycle is completed. 


3. Kinematic relations and mathematical model 


In the engine the spaces occupied by the working fluid are the 
hot volume, cold volume and regenerator volume. Hot volume is 
above the displacer top. Cold volume is between piston top and 
displacer bottom. Hot and cold volumes is calculated by using 
the kinematical relations. Fig. 2 shows the kinematic relations 
and schematic view of the engine with the rhombic-drive 
mechanism. 

The regenerator volume is the gap between the outer surface of 
the displacer and the inner surface of the cylinder and have con¬ 
stant value. The total volume of the engine was divided into 50 no¬ 
dal volumes: the hot volume (1), the regenerator volumes (48), and 
the cold volume (1). The following kinematic relations are used to 
calculate the values of cold and hot volumes: 



lrh 

2.666 



D • —1 

p r = sin 



R 


cr 


bh 


sin 9 




e c 


( R cr cos 9 - l rh cos p r + l r ) 


(ipr + R cr COS 9 + lrh COS p r + 




6h — U c ( Rcr COS 9 l r h COS /3 r -(- l^r T - ffrf) 

(4) 

Vc = ^c{Aq — A r ) 

(5) 

Vh = ChAh 

(6) 


All nodal volumes are open systems and consist of hot, cold and 
regenerator volumes. The working fluid friction losses between the 
nodal volumes are ignored and the following equation is used to 
calculate the pressure in the system: 


P 


m t R 


Y± _l_ Yl \ Yl \ I y_n 

T\ ' t 2 ' t 3 ' r T„ 



The temperature variation in nodal volumes is calculated using 
the first law of thermodynamics for open systems. The tempera¬ 
ture variations in nodal volumes are calculated by, 


AT _ [hjAj(T W j - Tj)At - AmjCJj + E, - pAVj] 

mjC p 



The following expression is used to calculate the enthalpy flow 
in or out of the nodal volumes: 


Ei = -C 


Ti + T i+ 1 


n 


P 


E Am J - c 

j ='+1 


T, i + T, . 

p —9 - z2 Am i 

j=i 



Hot and cold source temperatures are assumed to be constant. 
The surface temperature of the flow passage is supposed to be con¬ 
stant depending on time. The temperature changes linearly in the 
flow direction from cold source to hot source. Numerical solutions 
and detailed information of nodal analysis were given by Karabulut 
et al. [21]. 


4. Experimental setup 

The experimental setup is shown in Fig. 3. A Prony type dyna¬ 
mometer with an accuracy of 0.003 Nm was used to measure en- 
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gine load. The dynamometer consisted of a brake disk, a load cell 
and a torque arm. The engine speed was measured by a digital 
tachometer, ENDA ETS1410, with 1 rpm accuracy. The tempera¬ 
tures were measured by an infrared thermometer with ±2 °C accu¬ 
racy and -50 to 1600 °C measurement range. The experiments 
were carried out at a displacer cylinder hot-end temperature of 
673 K. Heat was supplied by an LPG burner. 


5. Results and discussion 

For the thermodynamic analysis, an engine speed of 500 rpm, 
cold source temperature of 300 K, and hot source temperature of 
773 K were determined. The inner heat transfer surface area of 
the displacer cylinder was enlarged by grooving 90 axial slots with 
1.5 mm width and 3 mm depth. For smooth and grooved cylinders, 
regenerator volumes were determined as 30 cm 3 and 90 cm 3 , and 
heat transfer surface areas were determined as 480 cm 2 and 
1200 cm 2 , respectively. 

Fig. 4 shows the cold volume, hot volume and total volume 
changes for grooved and smooth displacer cylinder. The variations 
of hot and cold volumes are the same for the grooved and smooth 
displacer cylinders. Channels inside the cylinder affect the total 
volume. The regenerator volume of the grooved cylinder is 
60 cm 3 higher compared to the smooth displacer cylinder. There¬ 
fore, the total volume of the engine with the grooved displacer cyl¬ 
inder is higher than that of the engine with the smooth displacer 
cylinder. As seen in Fig. 4, the interval of the crank angle from 
45° to 135° corresponds to the expansion process of the working 
fluid at constant temperature. The minimum cold volume appears 
at a crank angle of 90°. Since the cold volume remains at a mini¬ 
mum value in this crank angle range, the working fluid expands 
in the hot volume. The interval of crank angle from 135° to 225° 
carries out the cooling process - at constant volume. In this pro¬ 
cess, hot volume decreases and cold volume increases and the total 
volume remains constant. The compression process occurs at con¬ 
stant temperature in the 225-315° range of the crank angle. When 
the compression process starts at 225° crank angle, the cold and 
hot volumes are equal. During this process the hot volume de¬ 
creases. The cold volume increases up to the crank angle of 270°, 


-Total Volume (smooth ) -- Total Volume (grooved ) 

- Cold Volume Hot Volume 



Crank Angle (Degree) 

Fig. 4. Comparison of cold volume, hot volume and total volume variations 
obtained with grooved and smooth displacer cylinders. 

and then decreases. The heating process occurs in the range of 
315-45° crank angle at constant volume. 

p-V diagrams of the engine with grooved and smooth displacer 
cylinders at a charge pressure of 2 bar are given in Fig. 5 for 300 W/ 
m 2 K convective heat transfer coefficient and isothermal condi¬ 
tions. The dead volume of the engine with the grooved displacer 
cylinder engine is higher than the dead volume of the engine with 
the smooth displacer cylinder. The compression ratios of the en¬ 
gines with grooved and smooth displacer cylinders were calculated 
as 2.52 and 3.07, respectively. For a convective heat transfer coef¬ 
ficient of 300 W/m 2 K, engines with grooved and smooth displacer 
cylinders produce work of 21.37 and 18.4J, respectively. The cyclic 
top pressure of the engine with the grooved displacer cylinder is 
lower compared to the engine with the smooth displacer cylinder 
and its working fluid leakage is less than that of the smooth dis¬ 
placer cylinder. Due to the increase of the heat transfer area, the 
work obtained by the grooved displacer cylinder is higher com¬ 
pared to that of the smooth displacer cylinder. 

Fig. 6 shows variations of work with working fluid mass for con¬ 
vective heat transfer coefficients of 200, 300, and 400 W/m 2 K. 


LPG Heater 



Fig. 3. Schematic view of the experimental setup. 
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Fig. 5. Comparison of p-V diagrams obtained with grooved and smooth displacer 
cylinders. 


While the maximum values of work obtained for the smooth dis¬ 
placer cylinder with convective heat transfer coefficients of 200, 
300, and 400 W/m 2 K were 10.34, 15.52, and 20.63 J, respectively, 
the corresponding values obtained with the use of the grooved dis¬ 
placer cylinder were 14.09, 23.92, and 31.88 J. The engine work in¬ 
creases with increases in the working fluid mass. However, after a 
certain working fluid mass, the engine work decreases due to an 
insufficient heat transfer surface area. The engine work increases 
when using the grooved displacer cylinder due to its high heat 
transfer surface area. 

The variations of engine power with engine speed for different 
convective heat transfer coefficients are given in Fig. 7. Maximum 
engine powers for grooved and smooth displacer cylinders were 
obtained at an engine speed of 500 rpm. The maximum engine 
powers obtained with the smooth displacer cylinder were 86.30, 
129.48 and 172.63 W, respectively, for convective heat transfer 
coefficients of 200, 300, and 400 W/m 2 K. The maximum engine 
powers obtained with the grooved displacer cylinder were 
131.94, 199.40 and 265.87 W, respectively, for convective heat 
transfer coefficients of 200,300, and 400 W/m 2 K. Engine power in¬ 
creased up to an engine speed of 500 rpm and then decreased. The 
heat-exchange time decreased with increases in engine speed. 
Therefore, the amount of heat that is transferred to the working 
fluid becomes insufficient. 

The variations in engine power and torque obtained with the 
use of the smooth displacer cylinder are given in Fig. 8. The tests 
were conducted at a hot-end temperature of 673 I< and a cold- 


- 200 W/m 2 K (smooth ) 

- 300 W/m 2 K (smooth ) 

-400 W/m 2 K (smooth ) 


- 200 W/m 2 K(grooved) 

- 300 W/m 2 K (grooved) 

-400 W/m 2 K (grooved) 



Fig. 7. The variation of engine power with engine speed for different convective 
heat transfer coefficients. 
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a> 
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Fig. 8. The variation of engine power and torque with engine speed for the smooth 
displacer cylinder. 


end temperature of 300 K, with atmospheric air as the working 
fluid. Maximum engine power and cyclic work were found to be 
43.2 W and 6.35 J, respectively, at an engine speed of 300 rpm. 
The test conditions were carried out in the nodal analysis program 


-Isothermal - 200 W/m 2 K - 300 W/m 2 K 

- 400 W/m 2 K ▲ experimental 




Fig. 6. The variation of work per cycle with working fluid for different convective 
heat transfer coefficients. 


Fig. 9. Comparison of theoretical and experimental results obtained with atmo¬ 
spheric condition. 
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and the convective heat transfer coefficient was estimated as 
104 W/m 2 K. 

The comparison of experimental results and theoretical results 
is shown in Fig. 9. The values of cyclic work were calculated by iso¬ 
thermal and nodal analysis conducted with convective heat trans¬ 
fer coefficients of 200, 300, and 400 W/m 2 K. The difference 
between the experimental and theoretical results was quite small 
at low engine speeds. The difference between the experimental 
and theoretical values of cyclic work increased with engine speed. 
In the mathematical model, heat transfer between the working 
fluid and solid surfaces of the engine is taken into consideration. 
The heat-exchange time decreases with increasing engine speed. 
In addition, viscous flow losses and friction losses increase. These 
problems cause differences between the mathematical model and 
experimental results. 


6. Conclusions 

In this study, the effect of the grooved displacer cylinder on en¬ 
gine performance was investigated from kinematic and thermody¬ 
namic points of view. The maximum engine powers obtained with 
the use of the smooth displacer cylinder were 86.30, 129.48, and 
172.63 W, for convective heat transfer coefficients of 200, 300, 
and 400 W/m 2 K, respectively, while the corresponding engine 
powers obtained with the grooved displacer cylinder were 
131.94,199.40, and 265.87 W. Since both the heat transfer surface 
area and the regenerator volume increased with the use of the 
grooved displacer cylinder, engine performance was increased sig¬ 
nificantly. The maximum engine power obtained was 43.2 W at 
408 rpm in experiments conducted under atmospheric conditions 
using the smooth displacer cylinder. By nodal analysis conducted 
under experimental conditions, the convective heat transfer coeffi¬ 
cient was estimated to be 104 W/m 2 K. 
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